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Abstract 

We show that the factorization assumption in colour-suppressed B meson decays 
is not ruled out by experimental data on S — > K{K*) + J/^(^'). The problem 
previously pointed out might be due to an inadequate choice of hadronic form 
factors. 

Within the Isgur-Wise SU(2) heavy flavour symmetry framework, we search for 
possible gr^-dependences of form factors that satisfy both the large longitudinal 
polarization pL observed in B ^ K* + J/^ and the relatively small ratio of rates 
Rj/^ = T{B -^K* + J/^)/r{B + J/m). 

We find out that the puzzle could be essentially understood if the Ai{q^) form factor 
is frankly decreasing ( instead of being almost constant or increasing as commonly 

assumed ). 

Of course, the possibility of understanding experimental data is not necessarily a 
proof of factorization. 



PACS index 13. 25. Hw, 14. 40. Nd 
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I Introduction 

In a recent letter 0], Kamal and two of us (M.G. and X.Y.P.) have shown, within the factoriza- 
tion approach, the failure of commonly used B K{K*) form factors in explaining recent data 
on B —>■ J/'^ + K{K*) decays. The main problem is a simultaneous fit of the large fractional 
longitudinal polarization pi in B ^ J/'^ + K* decay and of the relatively small ratio of rates 
Rj/ifi between J/\E' + K* and J/\I' + i^" in the final states. We concluded that this difficulty in 
understanding experimental data might be due to a failure of the factorization method or to a 
wrong choice of hadronic form factors or both. 

Such an analysis has been independently performed by Aleksan, Le Yaouanc, Oliver, Pene 
and Raynal [0] who also found difficulties in explaining both and -Rj/-i., in spite of their large 
choice of heavy to light hadronic form factors consistent with asymptotic scaling law [Q . 

In our previous work [|l[|, in addition to our exploration of the usual B — > K{K*) form 
factors available in the literature, we also related the B K{K*) to the D K{K*) form 
factors using the SU(2) fiavour symmetry between the b and c quarks as first proposed by 
Isgur and Wise [Q. The input data are the hadronic form factors in the D sector at = 
as extracted from the semi-leptonic D — K{K*) + + ui decays. In such experiments, the 

distribution has not been measured and the analysis of experimental data has been made 
assuming monopole g^-dependence for all the D K(K*) form factors. For that reason in [0], 
we have also used monopole forms in the B sector. The resulting B K{K*) form factors 
obtained in this way are also unable to explain simultaneously experimental data on p^ and 

Our method, based on the Isgur- Wise relations, has been subsequently adopted by Cheng 
and Tseng [Q who considered various types of dependences for the hadronic form factors. 
However their model still encounters difficulties in reproducing correctly experimental data. 

The purpose of this paper is to make a purely phenomenological investigation of the possible 
g^-dependence — we shall call scenario — of the hadronic form factors in the B sector such that, 
assuming factorization and using the Isgur- Wise relations together with the latest data at 
g^ = in the D sector 0, we can obtain a good fit for both pL and -Rj/*. 

Some preliminary remarks are in order. We are aware of the fact that the values at g^ = 
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of the D K{K*) form factors have been extracted from semi-leptonic decay experiments 
assuming a monopole g^-dependence for all hadronic form factors. This ansatz is certainly 
inconsistent will all theoretical expectations coming, for instance, from QCD sum rules [^, from 
lattice gauge calculations as well as from scaling law of heavy flavours |^, || . A correct 
procedure would be to reanalyze the triple angular distribution fit in D semi-leptonic decay, 
with different scenarios, in order to evaluate the sensitivity to the scenarios of the normalization 
at = of the form factors. Such a study has not yet been done by experimentalists. Due to 
the limited range of in the D semi-leptonic decays, it was implicitly assumed that the values 
at = of the form factors could well be insensitive to the scenarios. To clarify and settle the 
issue, the cleanest information would come from measurements of the distributions for the 
rates and for the various polarizations in the semi-leptonic D sector. We are still far from such 
an ideal situation and for the time being, the only pragmatic way is to use the results quoted 
in with errors included for the values at g^ = of the form factors. 

We propose, in this paper, four types of scenarios for each of the B — > K{K*) hadronic 
form factors Fi,Ai,A2 and V in the Bauer, Stech and Wirbel (BSW henceforth) notation p. 

The g^ dependences, taken as (1 — g^/A^)~", are applied indiscriminately to all of these 
form factors. The algebraic integer n symbolically represents nF,ni,n2 and ny associated 
respectively to Fi,Ai,A2 and V. These integers n can take four values corresponding to four 
types of scenarios mentioned above : —1 for a linear dependence, for a constant, +1 for a 
monopole and +2 for a dipole. The respective pole masses Ap, Ai, A2 and Ay for Fi,Ai, A2 and 
V are treated as phenomenological parameters. Being related, in some way, to bound states 
of the bs system, we impose to these parameters the physical constraint to be in the range 
(5 — 6) GeV. Such a requirement is satisfied by the pole masses of the BSW model [§. 

We now summarize the results of our finding : 
The experimental data on pi and -Rj/-i. indeed can been fitted for three seenarios corresponding 
to three possibilities n2 = 2, 1, for ^* together with : 

i) ni = — 1 for a linear decreasing with g^ of Af^* 

ii) ny = +2 for a dipole increasing with g^ of V^^* 

iii) Up = +1 for a monopole increasing with g^ of 

pBK 

For a given selected scenario we have a non empty allowed domain in the A^, Ai, A2, Ay param- 
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eter space. Therefore we obtain hadronic form factors for B K{K*) reproducing correctly 
(within experimental errors) pi and -Rj/* with the parameters A^,Ai,A2,Ay physically ac- 
ceptable. We now easily understand why previous attempts fl], ||, Q were unsuccessful, mainly 
because the decrease with of the form factor Ai(g^) has never been seriously considered. Let 
us emphasize however that such an unusual behaviour has already been obtained by Narison 
in the QCD sum rule approach. Of course our result is not a proof of factorization in the B 
sector. It only makes wrong the statement that the difficulty in fitting simultaneously pi and 
Rj/ifi implies that factorization breaks down in colour-suppressed B decays. 

This paper is organized as follows. In part II we discuss in some detail the Isgur-Wise 
relations and, in particular, the consistency of scenarios in the B and D sectors as well as 
the relations between the parameters Aj {j = 1,2,V,F) in both sectors. The case of the form 
factors Fq and Aq, associated to the spin zero part of the currents is equally discussed. 

In part III we give the consequences of factorization for the decay amplitudes B — > K{K*) + 
{r]c, ■J/'^,'^') which are colour-suppressed processes. We study the kinematics and we review 
the available experimental data for these decay modes. 

The part IV is devoted to the decay modes B K(K*) + ^' , in which some scenario- 
independent results can be obtained. The left-right asymmetry A^^ji between the two transverse 
polarizations in B K* + ^' is found to be large and close to its maximal value. The fractional 
longitudinal polarization turns out to be a slowly varying function of A2 and the ratio of 
rates R^' a function of A2 and Ap. The result for i?^/ is consistent with experiment 
Our prediction for p'^ is compared with that of Kamal and Sanda who use seven different 
scenarios. 

The part V is the central part of this paper being related to the decay modes B K{K*) + 
J/^. The study of pl and Rj/^ allows us to select only three surviving scenarios among the 
4^ = 64 possible cases and to constraint the A^, Ai, A2, Ay parameter space. 

In part VI we make predictions for the decay modes B K{K*) + rjc for which, unfortu- 
nately no experimental data are available. 

The comparison of various charmonium states involves different leptonic decay constants 
/<^') fj/^ and fr,^. This problem is studied in part VII in relation with previous works and 
|[10| and with experimental data available only for the ratios of ^' and J /'^ production. 
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The B K* vector and axial vector form factors studied here can be related under reason- 
able assumptions to the tensor and pseudotensor B —>■ K* form factors describing the radiative 
decay B — > K* + 7 recently observed. This problem is briefly discussed in part VIII. 

Finally in part IX, we come back to the D sector in the light of results obtained in the B 
sector. Of course, the hadronic form factors Fi, Ai, ^2, follow the same scenarios in the B 
and D sectors with poles masses related by Equations ( |II.24]) and (p.37|) . We determine the 



normalized distributions for semi-leptonic decays D — > Ki~^iyi, D — > K* l^vi and for this last 
mode the integrated longitudinal polarization pf^ and the left-right asymmetry Af^p^ between 
transverse polarizations {si stands for semi-leptonic). We also study the hadronic two body 
decay modes D° — > K^{K^*) + 7r+(p+) which are colour favoured and we show that simple 
factorization assumption fails in the D sector where, in addition, final state strong interactions 
play an important role. 

Discussions and critical remarks on the results are given in the conclusion. 

II The Isgur-Wise relations due to SU(2) heavy flavour 
symmetry 

1°) The SU(2) flavour symmetry between the heavy b and c quarks allows us to derive relations 
between the B K{K*) and D —>■ K{K*) hadronic form factors at the same velocity transfer 
though at different momentum transfers. Calling ts [to) the value of the squared momentum 
transfer for B{D) form factors, we obtain the following kinematic relations : 

Vb-VK = Vc-VK or met B -rnbtD - {rrib- rric) {rrib rric- m]^) = (H-l) 
Vb ■ Vk* = Vc ■ Vk* or rrict*^ — rribfj^ ~ {i^b — rric) (nib mc — rn\,) = (II-2) 

In particular, at zero squared momentum transfer in the D sector, t^i = = 0, we get 

l){mb rUc - m]^) (II. 3) 



tB 


— fO — 

— Ib — 


^rrib 




rric 


f* 


— Ib — 


rUb 
rric 



40 (11.4) 



and useful relations such as 



771 

tD = — (tB - t^) (II.5) 

rrib 
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t 



D 



— it*B 

nib 



t 



B . 



(11.6) 



The knowledge of the hadronic form factors at = in the D sector will determine the 
hadronic form factors in the B sector at = for the B ^ K case and at = t*^ for the 
B ^ K* case. 

The values of and t*^ depend on the quark masses and rric- We choose, in this paper, 
rrib = 4.7 GeV, rric = 1.45 GeV and get : 

14.7287 GeV^ 



f° 



A.*0 

'^B 



13.4933 GeV'^ 



(11.7) 
(11.8) 



2°) We first consider the case of B ^ K and D ^ K form factors. The matrix elements of 
the weak current involve two form factors /+ and /_ defined by 

<K\J,\D> = (pd+Pk), f^''iq') + iPD-pK), f^^^iq') (II.9) 
<K\J^\B> = (p^+p^)^ f^^{q') + {p^-p^)^ f^^{q') (H.IO) 



where q = Pb,d - Pk- 

The Isgur-Wise relations are written as 



Che 
Cbc 



nib 



Irrib 



iU - /-)^^(tD) 



nir. 



or equivalentlly 



be 



mb + rric oKr. ^ 



(11.11) 

(11.12) 
(11.13) 



_2y/m^ 2y/rfiymc 
where Cbc = [^^{^ ]"^/^^ is the QCD correction factor. | 

In the BSW basis , the spin one and the spin zero parts of the weak current are separated 
and two new form factors Fi and Fq are defined 



Friq'' 



iq' 



f^iq') 

/r(^^) + 



m. 



m 



f-^'iq') 



(11.14) 
(11.15) 



K 



^In this paper, we use C^c = ] ^^^^ — 1.135 which has been obtained by using the recent world 



averaged values for A-p-j . —^^^^ 
rrib = 4.7 GeV, = 1.45 GeV. 



= (225 ± 85) MeV, Ai^ = (325 ± 110) MeV, and for the quark masses : 
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where P = B or D. 

Let us define the ratio of form factors /_ and by 



^^ir) = --7pwrJ^ P = B,D (11.16) 



Using Eq.( |11.13| ), we obtain a relation between the two quantities fJ^^its) and fi^{t£)) 



1 + a fi^ito) 

where the constant a depends only on the quark masses 



l^'its) = '^^"t\ (11.17) 



a = ^^^t^ (11.18) 

rUb + rric 



The spin one function 

pBK 

is now related to f'^^ and fi^, and we obtain from Eq.( [11.13D 



Frits) = ^^^^^ [l + a /i^(tz.)] Fnto) (11.19) 
and for the spin zero function F(^^ , we get from Eq.( [11.15| ) 

(g^) = [1 - /.^(g^)] Fne) (11.20) 

3°) As explained previously, we shall use in the D sector the values of the hadronic form 
factors at = as coming from semi-leptonic experimental data. Because of the normalization 
constraint F[^^{0) = F^^{0) only f^lO) is known and we cannot have, in that way, any 
information on f^^{0) and /i^(0) f\ . 

We now make an assumption which is natural in the framework of the SU(2) heavy flavour 
symmetry. If f^iq"^) and f^^{q^) have the same type of g^-dependence, no matter how it 
is, then the ratio is constant and, using the Isgur-Wise relations ( [11.11| ) and ( [11.12| ) we 
easily see that the same property extends to the B sector and , in particular, the ratio is a 
constant related to by Eq.( [lI.17D . Of course, Ff^lq"^) and F^^^lq"^) have the same type of 
g^ dependence. 

For instance, if F['^{q^) is written in the form 

pDK^2^ = ^P""^^) (11.21) 

^ A^;] 

§ In principle, f?^ {0} and /i'^(O) could be measured in D ^ K + + i/^ namely by looking at polarized 
muon. 
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where np is some algebraic integer, then, using Eq.( [il.l9D we obtain a similar expression for 
with the same np. 

Furthermore a relation between the pole masses Kp and h^oF has been previously given in 
Ref. and it comes from the identity 



+ (".23) 

The result is 



rric A| - mfc A|,^ = = (m^ - mc){mh rUc - m^) (11.24) 

The values at = of the form factor Fi in the B and D sectors are also related by Eq.( |ll.l9| ). 



The result depends on np and is given by 

4°) An interesting scenario for suggested by many theoretical studies 0, |ll|] as well 

as supported by experimental data [^, |12[ is a monopole dependence n^? = 1 with a pole mass 



Adf in the 2 GeV^ region. While F^^^ and Ff^ increase with by a monopole type, from 
Eg. ( [II. 20]) we easily see that the form factors F^^ and F^^ have a different behaviour due 
to a supplementary linearly decreasing factor in Eq.( |11.2Cl| ) 

[1-Z:t^ for Fo^^ 



m R — m 



[1- , /i^] for Fo^^ 

In fact, for a particular value of fi^ or fi^, this factor can exactly cancel the pole of Fi in the 
B or in the D sector, thus making Fq constant. This is the case when 

2 2 2 2 

B '^B ~ D ~ /TT r,«^ 

/i = -2 or /i = -2 (11.26) 

^^F ^^DF 

In this way, the value of fi^{fi^) is related to the pole mass Ap{Ap)p). However such a situation, 
in general, does not occur in both B and D sectors because A|. and A|,^ are not independent 
from each other due to Eq. ( [11. 24] ). 
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If we impose the constraints ( |11.26| ) for fi^ and fi^, we can determine the pole masses A^. 
and A'jjp. The relations ([1I.17|) , (|II.24| ) and ( [1I.26|) lead to a second order equation for Aj^p, the 
roots of which are real. However one root ( negative) is physically unacceptable, and the other 
one ( positive ) gives A^p within the (2 — 3) GeV range, corresponding to Ai? within the (5 — 

6) GeV range and the value of fi^ is in the neighbourhood of 1. 

The calculation involves the b and c quark masses : with = 4.7 GeV and rric = 1.45 
GeV, the hadronic form factor Fq is constant in both B and D sectors for the following pole 
mass values A^f = 2.32 GeV, Ap = 5.67 GeV and for fi^ and fi^, we obtain fi^ = 0.8585 and 

= 0.6041. 

5°) We do not have, to our knowledge, a theoretical estimate of the form factors f^^{0) and 
f?^{0). The heavy quark symmetry limit can not be applied for heavy quark to light quark, 
b —>■ s, transitions. Therefore we shall use, in the B sector, a model suggested by theoretical 
considerations 0, 13 1 where is independent of and Ff^ has a monopole dependence 
with a pole mass Ap. The parameter is then a function of A'j^ as given in Eq. (p.26|) and 
the parameter /i^ is known from fi^ by inverting Eq.( [ll.l'^ ) 

where a = 0.5285 for our choice of quark masses. 

We plot in Figure 1 the ratios fi^ and fi^ as funtions of A|. for values of A^t^ in the (4 — 

7) GeV range. They are parts of hyperbolae with asymptotia parallel to the A|. and yU axes. 
We plot, by a straight line, Eq. (|II.24| ) which relates the pole masses in the B and D sectors. 
On the same figure, we also show, with dotted points, the quantity Ji^ which corresponds to a 
constant Fq form factor in the D sector. 

6°) We now study the case of the B —>■ K* and D K* form factors. The matrix elements 
of the weak current involve four form factors, /, g, a+ and a_ in the Isgur-Wise basis, or Ai, V, A2 
and Aq in the BSW basis. 

The Isgur-Wise relations for / and g are very simple P] 



f^^\f^) = ac /^^*(t^,) (11.28) 



g'^^^fB) = g^^\tl) (11.29) 
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Using the relations between the Isgur-Wise and BSW bases : 

1 



/ 

we get 



PK* 



PK* 



{mp + rriK*) 
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PK* 



{mp + rriK*] 



V 



PK* 



be 



A^^\ep) 



TTic mp + niK* 



P = B,D (11.30) 



(11.31) 



V mf, mp + ttlk' 



In particular, defining the ratio 



Ai{q' 



(11.32) 



(11.33) 



the values of y in the B and D sectors, and y^ are related by : 

"^6 \mp + mK*J 

The dependence with respect to of Ai and \^ in the B and D sectors is preserved by the 
Isgur-Wise relations (p.31|) and ( [1I.32|) . We choose the form factors A^^* {q^) and V^^' {q^) 
such as : 

Af^*(0) 



V^^^*(g2) 



[1 - S]- 



(0) 

[1 - I;]- 



(11.35) 



(11.36) 



and analogous expressions in the D sectors with pole masses A^i and Kpv related to Ai and 
Ay by a formula similar to Eg. ( [II. 24]) 



rric A| 



rrib Ap = rric tp = {nib - mc){mb rric - m\,) 



(11.37) 



where A^ = Ai or Ay and Kp = Api or Apy. The algebraic integers rii and ny will be 
restricted to the values —1,0, 1 and 2. The normalizations at = of these form factors in 
the B and D sectors are related using Eqs. (|II.31|) and ( 11.32 ) 

Af^'*(0) 



Af^*(0) 



Chr \ 



nib mp + rriK* 



rric + rriK* 



J-*0 

^ AV 



(11.38) 
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7°) For the two other form factors a+ and a_ describing the P K* transition, the situation 
is formally similar to and /_ in the P K case previously considered. The Isgur-Wise 



relations are written as 



(a+ + a„)^^*(t^) =C,A — )^ (a+ + a_)^^*(tj,) (11.40) 
(a+-a_)^^*(t^) = ac {—\" (a+-a_)^^*(tB) (11.41) 



\ 3/2 



from which we deduce 

1 „ (rrir. 



2 "'^ \m. 



rrib 



;i + ^)af^*(tB)-(l-^)aJ^*(tB) 



(11.42) 



In the BSW basis p], the spin one and the spin zero parts of the weak current are separated, 
with Ai, V and A2 for the spin one part and Aq for the spin zero part. The form factors A2 
and Ao are linear combinations of the form factors a+, a_ in the Isgur-Wise basis : 

= -{mp + niK*) a^'^'iq') (11.43) 

APK*^^2^ ^ mp + rriK* j^pk* ^^2^ _ mp- rriK* ^pk* ^^2^ ^ a^^*(g^) (11.44) 

2 Tfix* 2 TTix* 2 nix* 

where P = B or D. 

At = the normalization of Aq is determined by Ai(0) and ^2(0) : 

(0) = ^^^^^^ Af-(o) - ^^^^^^ (0), p = 5,/^ (11.45) 

and we can not have any information on the value at = of the form factor a_. |^ Let us 
now define the ratio of the form factors a_ and a+ by : 



^''(^') = - ^^pF^ P = B,D (11.46) 



Using Eq.(|L42) we obtain a relation between the two quantities X^{t*^) and A'^(t^) similar 
to Eq. (pT7D 



^similar to the f-{q^) case, the form factor a-{q^) can be measured m D ^ if + pt^ 
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with o defined in Eq. (|11.18| ). 

The function Af^*(t^) is now expressed in terms of A^^* {t*j^ and A'^(t^ 
Using Eqs.(|L42D and ( [11431 ), we get 



B) 



rrir 



mB + rriK* r. , ^ \D(+*\^ aDK*(+* 



nib/ rriD + mK* 



D) 



(11.48) 



With Eqs.(|L4|) and (|L4D we obtain the spin zero function 



BK* 



2 

Defining the ratio : 



2 ttlk* 



[1 



ml 



m 



K* 



Ml 

AAq^ 



A^(g^)] A^^ {q') (11.49) 



(11.50) 



and using Eqs.( ii.3"T|) and ( if. 48 ), the values of x in the B and D sectors, and x^ , are related 
by 



1 m, 



2 mfe 



f^B + rriK* 



+ a A^(tl,)] x^{fj,) 



(11.51) 



8") We now make an assumption for a+ and a_ similar to the one made previously for /+ 
and /_. If a^^* {q^) and a^^* {q^) have the same type of g^-dependence, no matter how it is, 
then the ratio A^ is constant and using the Isgur-Wise relations ([11.40| ) and ([11.41| ), the same 
property extends to the B sector with, in particular, a constant value for the ratio A"^ related to 
A^ by Eq. ( |II.47]) . Of course, A^^* {q^) and A!^^* {q^) have also the same type of g^-dependence. 

The form factor A^'^*(g^) is written in the form 

XO) 



Ar{q' 



n2 



(11.52) 



and the normalizations at = of A2 in the B and D sectors are related using Eq.( [11.48D 



(0) 



(0) 



rrib 



be 



1 + — 

nib 



m-B + niK* 
niD + niK* 



[l + a\ 



j.*o 

^ AV 



(11.53) 



The second term in Eq.( [11.49[ ) for Aq is written as 



^B + mx* 
2 nriK* 



A?^^{q' 



rriB - niK* 
2 niK* 



[1 



m 



B 



ni 



A«] 



Af^*(0) 



K* 



"2 



(11.54) 



and it exhibits a different q^ behaviour compared to A2 due to a supplementary linearly de- 



creasing factor (1 



). For a particular value of A^ or A^, this factor cancels exactly 



D 
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one power of A^'^ or A^^*. This is the case when 



m-R — ml', , n m'^n — rn? 



A2 A2 



In this way the value of A'^(A'^) is related to the pole mass A2(A£)2)- However such a 
situation does not occur, in general, in both B and D sectors because and A|,2 are not 



independent from each other due to a relation similar to Eq.( II.37 ). Conversely if we impose 
the two constraints ( |II.55|) , we can determine the pole masses A^ and A|)2. The relations ( [1I.37| ) 
and ( [11.471 ) lead to a second order equation for A|,2- However one root is negative and must be 
rejected. The second one gives A£,2 in the (2 — 3) GeV range corresponding to A2 in the (5 — 

6) GeV range. The value of is not far from 1, the heavy quark symmetry prediction. Using 
as previously rjif, = 4.7 GeV and rric = 1.45 GeV, we obtain Ad2 = 2.59 GeV, A2 = 5.93 GeV, 
XB = 0.7698 and A^ = 0.4068. 

9°) In the absence of a theoretical estimate for a^^* (0) or a^^*{0), we shall use in the B 
sector a model where A^ is related to the pole mass A2 by Eq. ([II.55|) . The parameter A"^ is 
then computed from A^ by inverting Eq.( |II.47|) 

- 1& 

We plot in Figure 2 the ratios A"^ and A^ as functions of A2, for values of A2 in the range (4- 

7) GeV. They are parts of hyperbolae with asymptotia parallel to the A2 and A axes. Eq.( |11.37| ) 
which relates the pole masses in the B and D sectors is represented by a straight line in Figure 
2. On the same Figure we also show, with dotted points, the quantity A^ which corresponds 
to the second relation ( [1I.55| ). 



Ill Factorization and Kinematics, Experimental Data 

1°) The two-body decays of the charged and neutral B mesons discussed in this paper are 
described, at the tree level, by the colour-suppressed diagram of Figure 3. Of course penguin 
diagrams also contribute to these decays at the one loop level. However the colourless char- 
monium states cc have to be excited from the vacuum and for which two or three gluons are 
needed. For that reason the penguin contribution will be neglected in this paper. 
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2°) We consider the decay modes 

(III.l) 

K°{K*°) +r],{J/'^, 

and we compute the decay amphtudes assuming factorization. We obtain an expression of the 
form 

<cc + sq\T\bq> oc < cc\J^\0 > < sq\J^\B > (III.2) 

The first factor in the right hand side of Eq.(|TL^) involves the leptonic decay constants frj^, fj/q, 
and for rjc, J/\l/ and respectively. The second factor is governed by the hadronic from 
factors for the B ^ K ot B K* transitions. As a consequence, the branching ratios have 
the following structure 

BR = BRo ■ I — ] ■ PS ■ FF (III.3) 

The scale BRq contains the Fermi coupling constant Gp, the Cabibbo-Kobayashi-Maskawa 
( CKM ) factors V^^Ks as indicated on Figure 3, the B meson life time tb and the BSW 
phenomenological constant 02 describing colour-suppressed processes 



BRn 



\Vcs\' ^ «2 f (111.4) 



V2 

Being interested only in ratios of decay widths, we shall not numerically compute BRq. 

The quantity PS is a dimensionless phase space factor depending only on masses of the 
involved particles. Because of the small K~^{K*~^) — K°[K*°) mass differences, the numerical 
value of PS is slightly different for B^ and B° decays. However the differences in these quan- 
tities are typically (9(10~^), hence we ignore the mass difference between charged and neutral 
strange mesons. We give results for B^ decay. 

The last factor FF depends on the hadronic form factors and it contains the dynamics of 
the weak decays. The results are : 

(a) B+^K+ + 7], : = 0.3265 , FF = \F^^ {mlJ\'^ (III.5) 

(6) 5+^i^+ + J/^ : P5 = 0.1296 , FF =\Ff^{m]/^)\^ (III.6) 
(c) B+^K+ + '$' : P5 = 0.0575 , PP = |Pf ^(m^,)P (HI.T) 
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(d) B+^K*++r], : P5 = 0.1218 , FF = |A^^*(m;j|2 (III.8) 

(e) 5+ ^ K*^ + J/^ : PS = 0.1399 , 

FF = |Af^'(m2,/^)p [(a - + 2 (1 + cV)] (ni.9) 
(/) fi+ ^ + ^' : P5 = 0.1407 , 

FF = |Af^*KOr [{a - fe'x? + 2 (1 + c'V')] (III.IO) 



The analytic expressions for a, 6, c are previously given in Ref. ||l[ and a ,b' , c are obtained 
respectively from a, b, c by the simple substitution m^> to mj/^. We get numerically : 

a = 3.1652 6= 1.3081 c = 0.4356 (III.ll) 

a = 2.0514 b' = 0.5538 c = 0.3092 (III.12) 



The ratios of form factors x, y, x' , y are defined by : 

x^x (m^/^) = -g^^— , y = y {mji^) = —^^—^ (III.13) 



In the i^'* + J/\I^ and i^* + modes, we have three possible polarization states, one is 
longitudinal and two are transverse for both final particles. We shall define two interesting 
quantities, the fractional longitudinal polarization : 

T{B K* + J/^) ' T{B ^ K* + ^') ^ ' 

and the left-right asymmetry : 

T{B ^ K* + J/^)— - T{B ^ K* + J/^)++ 



A 



LR 



T{B K* + J/-^)— + V{B ^K* + J/^)_ 
T{B ^K* + ^')__ - T{B -^K* + ^')++ 



(III.16) 



We get 



(a-te)2 + 2[l + c2y2] ' Pi^ ^a' -b'x'y + 2[l + c"y"] ^ ^ 
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A 



LR 



Icy 



1 + ^2 

We shall also define two ratios of rates 



X 



Icy 

l + c y 



R 



J/* 



T{B ^ K* + J/^) 



R 



3°) The experimental data for decay rates as averaged by PDG |^ are given in Table 1. 



(111.18) 



(111.19) 



Modes 


B+ 


B° 


ir+ ^' 


(0.69 ±0.31) • 10-3 


< 0.8 • 10-3 


K* + ^' 


< 3.0 ■ 10-3 


(1.4 ±0.9) ■ 10-3 




(1.02 ±0.14) ■ 10-3 


(0.75 ±0.21) ■ 10-3 


/r* + j/^ 


(1.7 ±0.5) ■ 10-3 


(1.58 ±0.28) ■ 10-3 



Table 1. 



The modes Krjc and K*r]c have not yet been observed experimentally. 
The ratios R^' and -Rj/^ can be estimated from the data of Table 1 and the results are 
shown in Table 2. 



Ratios 


B+ 


B" 


B^, B° combined 


R^' 


< 4.35 ±1.95 


> 1.75 ± 1.12 


2.03 ± 1.59 


Rj/^ 


1.67 ±0.54 


2.11 ±0.70 


1.83 ±0.43 



Table 2. 



A direct measurement of Rj/q, by CLEO 11 [|14| is consistent with our estimate given in the 
last column of Table 2, 

CLEO 11 = 1.71 ±0.34 

In what follows we shall use the constraint -Rj/* < 2.2. 

A second type of useful experimental data which turns out to be crucial is the fractional 
longitudinal polarization for the mode B —>■ K* ± J/\l' measured by 3 groups : 
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CLEO II H PL = 0.80 ± 0.08 ± 0.05 

CDF H pL = 0M±0.10t [J:i° 



ARGUS ig PL = 0.97 ±0.16 ±0.15 



Averaging these results by the standard weighted least-squares procedure, we obtain 

= 0.780 ± 0.073 

In what follows we shall use the one standard deviation lower limit pi > 0.7. 



IV The decay modes B K{K*) + ^' 

1°) We first consider the decay B K* + \E'' which is described by the three hadronic form 
factors Af^* , A^^* and V^^* taken at = m^, . In part II we have explained how the Isgur- 
Wise relations due to SU(2) heavy flavour symmetry allow one to compute the form factors in 
the B sector at = t*£ from the values of the form factors in the D sector at = 0. It turns 
out that the numerical value of t*£ = 13.4933 GeV^, as given by Eg. (|II.8|) , is remarkably close 
to m^/ = 13.5866 GeV^. Different plausible choices of rrib, rric satisfying mj, — rric = 3.4 ± 0.2 
GeV (as dictated by HQET scheme) yield similar results, t*£ is close to m'^t. It is then justified 
to neglect the variation of the form factors between t*£ and m^/ and we obtain from Eqs. (|II.31|) 
and dig) : 

^f^'K') = 0.8056 Af^*(0) (IV.l) 
\/^^'(m^,) = 1.2413 Cbc V'^^'iO) (IV.2) 

With the values in the D sector as given by the Particle Data Group ( PDG henceforth ) : 

Af^*(0) = 0.56 ±0.04 , l/^-^*(0) = 1.1±0.2 (IV.3) 

we get : 

^f^*(^l') = (0.4511 ±0.0322) Cbc (IV.4) 
V^^'iml,) = (1.3654 ±0.2483) Cbc (IV.5) 

In an analogous way, we have : 

y^'iml,) = 1.5409 7/^(0) (IV.6) 
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and with the PDG values §, y^{0) = 1.89 ± 0.25, we get 

y' = y^(m^,) = 2.9123 ±0.3852 (IV.7) 

The values of Af^ and V^^* at = m'^, are simply given by the Isgur-Wise relations from 
the values of A^^* and V^^ at = and the relative error is the same in both B and D 
sectors. 

For ^* and the situation is not so simple because of the presence of the parameter 

= [0.8121 + 0.4291 A^] Cbc Af^*(0) (1V.8) 
x = x^{ml,) = [1.0081 + 0.5328 A^] x^(0) (1V.9) 

Using the model introduced in Part II, A^ and A^ are functions of the pole mass A2 of the 
form factor Af^* and represented on Figure 2. Using the PDG values [^ : 

Af^*(0) = 0.40 ±0.08 , a;^(0) = 0.73 ±0.15 (IV.IO) 

we obtain the quantities Af'^*(m^,) and x^^ni'^,) as two functions of A2. They are represented 
on Figure 4 for A2 in the range (5 - 6) GeV. 

2°) The knowledge of y' = y^^m^,) determines the left-right asymmetry A'^^j^ previously 
defined. The result 

= 0.9945 ±0.0137 (IV.ll) 

shows that the dominant transverse amplitude has the helicity A = — 1. In the one standard 
deviation limit, A^^ji > 0.98. 

As a second consequence of the knowledge of y' , we can derive an upper bound for the 
longitudinal polarization : 

,2 

Pl — f2 1 on I /2 /2-| (IV. 12^ 

a ± 2[1 ± c y \ 
In the one standard deviation limit for y\ we obtain 

p'l < 0.5664 (IV.13) 

which is significantly smaller than the theoretical upper bound Pl < a /{a ± 2) = 0.678. 
These two results ( [1V.11| ) and ( |1V.13| ) are clearly scenario independent and they are direct 
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consequences of the Isgur-Wise SU(2) heavy flavour symmetry assuming the numerical value 
of y^{0) as given by experiment to be correct. 

3°) Furthermore not only the upper bound Eq.( LV.l"3 ) but also exact values for can be 



computed from x = x^{m'^,) and y' = y^{m?^i). It is a function of A2 represented on Figure 
5 for A2 in the (5 — 6) GeV range. The error on p^, Ap^, combines, in a quadratic way, the 
error on x and y' due to those on x^(0) and y^{0). We observe that p^ is a slowly increasing 
function of A2 which takes the value p'^ = 0.403 ± 0.042 for A2 = 6 GeV. In Figure 5 we also 
represent our one standard deviation upper bound ( [IV.13|) . 

Estimates for p^ have been obtained by Kamal and Santra ||^. However their method 
and results are different from ours. Seven scenarios are considered for relating the J/"^ and 
^' modes and the allowed domains in the x\y' plane for each scenario are limited by the 
constraint pl > 0.68. The upper bound for p^ found in ^ is the theoretical upper bound 
p'l max = a'^/(a'^ + 2) = 0.678. Their lower bound is slightly scenario-dependent and it varies 
from 0.48 to 0.55. In our case, the allowed domain for p^, with one standard deviation, is 
smaller than their results and we have : 0.286 < p^ < 0.445. 

Let us emphasize again that our prediction is scenario- independent. It is determined by 
using both the Isgur-Wise relation and the experimental values of a;'^(0), y^{0). It is interesting 
to compare the experimentally observed pi with the theoretical prediction for p^ : 

PL = 0.78 ± 0.073 , 0.286 < p'l < 0.445 

4°) Let us consider now the decay mode B K + "^' described by the hadronic form factor 
Ff^(m^,). The Isgur-Wise relation (|II.19|) gives Ff^{t%) in terms of F/^^(0) and Ap via the 
parameter p^ 

Ff^(t^) = [1.1779 + 0.6225 p^] Cb, F^^{0) (IV.14) 

We use the PDG value § ^(0) = 0.75 ± 0.03. 

In order to obtain Ff^{m'^i), we must extrapolate Ff^ from t% = 14.7288 GeV^ to m^, = 
13.5866. For that purpose we use, as explained in Part II, a monopole g^-dependence with a 
pole mass Ap. The result is shown on Figure 6 with Ap in the (5 — 6) GeV range. We observe 
that Ff^^m^,) is a slowly decreasing function of Ap. 
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5") We finally study the ratio of rates 

T{B ^ K* + <if' 



r{B ^ K + -^') 

Using the phase space factors given in Part III, we obtain 



/2 /2-, 



= 2.4465 (-'-tyf + all + c y . 

z 

where the ratio of form factors, z\ is defined by 

' KO- ^k.^) (IV.16) 

The ratio z\ like Ff^{m'^,), is a function of A|, and the ratio x is a function of A|. It follows 
that i?^' depends on both parameters Ap and A|. Restricting Ap and A2 in the range (5 — 6) 
GeV, we find the ratio R^' between 1.44 ± 0.28 (for A^ = A2 = 5 GeV) and 2.92 ± 0.54 (for 

Ap = A2 = 6 GeV). 

Both extreme values are obviously compatible with the experimental estimate given in Part 
III : 

(i?*')exp = 2.03 ±1.59 (IV.17) 

The ratio of rates R^' is represented on Figure 7 as a function of A2 and A^ both in the range 
(5-6) GeV. 



V The decay modes B K{K*) + J/^ 

1") We consider the decay mode B K*-\-J/'^ described by the three form factors Af^* {rrij^^), 
Af^*(mj/^) and V^^* (m'^j^^). These quantities are related to the values at t*^ of the same 
form factors , already obtained in Part IV, via the scenario dependent parameters cti, a2 and 
p. 

^r*K*) = «2^?^*(iS') (v.i) 
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and for their ratios 

X ^ a;«K/^) = px^ra (V.2) 

y ^ y^'im]/^) = q y^f^) (V.3) 

where 

p = — , q = — (V.4) 

We shall consider the scenarios qualitatively described in the Part II. For that purpose, we 
introduce the function r(A) defined by : 

A 2 j.*o 

KA) = j^^^ (V.5) 

and we get : 

ai = [r(Ai)r , a2 = HA2)r , /? = [r(Ayr- (V.6) 

For the algebraic integers ni,n2,ny, we shall consider the four cases rii = — l,0,+l,+2 with 
i = 1,2,V. On physical grounds, we impose to the pole masses Ai, A2, Ay to be inside a cube 
5 GeV <Ai<6 GeV. 

2°) We first study the scenario constraints due to the longitudinal polarization fraction pl in 
B —>■ K* + J/\I^. For the 4'^ = 64 possible triplets [rii, n2, ny], we compute pi using the values of 



X 



^{t*^) and y^{t*£) obtained in Eqs.(|V]|) and [W^ and restricting the pole mass parameters 



Ai, A2 and Ay inside the cube 5 GeV < Aj < 6 GeV . We impose the experimental constraint 
in the form + Ap^ > 0.7 where the error Ap/, is computed in quadrature [See Appendix 
A-1] from the errors on x and y due to the experimental errors of x^(0) and y^{0). After 
numerically scaning over the whole Ui and Aj spaces, finally our results can be summarized in 
the following : 

1. No solution is obtained when rii = 0, 1, 2 for the 16 values of the couple [n2, ny]- 

2. Solutions exist only when rii = —1, i.e., when the form factor Ai exhibits a linear decrease 
with g^. Of course, in this case, Ai is no more a pole mass but simply a slope coefficient 
and it is reasonable now to relax the constraint Ai < 6 GeV and to use only Ai > 5 GeV 
in order to exclude a too fast variation with of Ai. 
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3. We obtain solutions for only 4 triplets [ni,n2, ny] '■ 

[-1,2,2] ; [-1,1,2] ; [-1,0,2] ; [-1,2,1] 

The allowed domains in the Ai, A2, and Ay space are respectively represented on Figures 
8, 9, 10 and 11. 

4. In the four surviving triplets mentioned above, the maximal value of pl occurs at Ai = 
5 GeV, A2 6 GeV, Ay = 5 GeV and in the most favorable situation of two dipole q'^ 

dependence for A2 and V, we obtain = 0.7162 ± 0.0236. Therefore pi = 0.74 is the 
maximal value, within one standard deviation, we can get in our approach, considering 
only the quantity p^. 

5. It is interesting to notice that in the case of two monopolc g^-dcpcndence for A2 and V, 
the maximal value of pl is obtained at the point Ai = Ay = 5 GeV, A2 = 6 GeV with 
the result pl — 0.6635 ± 0.0339, e.g., a one standard deviation value very close to 0.7. 
However we consider this case as only marginal, since, in any way, it will be eliminated 
when the second quantity Rj/^ enters in the fit. 

3") The second quantity is the ratio of rates which has the form : 

. 1.0793 (a-6x)- + 2(l + c^,^) ^^^^ 

z 

where the ratio z is defined by : 

z = z (rrij/^) = -^j^^—^ (V.8) 

We impose the experimental constraint under the form — ARj/-q, < 2.2 where the theoret- 
ical error ARj/q, is computed, in quadrature [See Appendix A-3], from the experimental errors 
on x^{0), y^{Q) and ^-^(0). 

Our constraint pl + Api > 0.7 has selected four scenarios and the allowed domain at fixed 
Ay, A2 is defined by : 

Ai < Ai,MAx(Ay,A2) (V.9) 

These domains have been represented on Figures 8 to 11 for 5 GeV < Ay, A2 < 6 GeV. 
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On the other hand, the constraint on Rj/m imphes a lower hmit for Ai 

Ai > Ai,M/iv(Ay,A2,Air) (V.IO) 



Physical values of Ai exist when and only when the lower limit ( [V.10|) is smaller than the upper 
hmit (|VJ| ). 

At fixed Ay, A2, the quantity Ai. min is an increasing function of Ai? and we shall restrict 
Kp ^ 5 GeV . Therefore at fixed Ay, A2, the physical domain for Ai is defined by 

Ai, M/7v(Av, A2, Ai. = 5 GeV) < Ai < Ai, MAxi^v, A2) (V.ll) 

After a numerical scanning, we do not obtain solution satifying the constraint Rj/is, — 
Ai?j/^ < 2.2 with the scenario [ni,n2,'nv] = [—1,2,1]. For the three remaining scenarios 
[ni,n2,nv] = [— 1, 2, 2], [— 1, 1, 2] and [—1,0,2], the physical regions are then represented on 
Figures 12, 13 and 14. 

At fixed Av,A2, we have 

5 GeV < Af{Av,A2) < Ap, MAx{Av,A2) (V.12) 

where Ap, max is determined by 

Ai^ min[Av, A2, Ap^ MAxiAv, A2)] = Ai,MAx(Av,A2) (V.13) 

The quantity Ap^ max{Av, A2) is represented on Figures 15, 16 and 17 for the three surviving 
scenarios. 

4°) Comment and Illustration 

Starting with 64 scenarios for the g^-dependence of the hadronic form factors Ai,A2,V, 
[ni,n2, ny] with Ui = —1, 0, 1, 2, only survive 3 scenarios [—1, n2, 2] with n2 = 0, 1, 2 for which 
it is possible to find a non empty domain in the parameter space Ai, A2, Av,Ap such that both 
experimental constraints pi + Ap^ > 0.7 and -Rj/* — ARj/^ < 2.2 are simultaneously satisfied. 
The hadronic form factors F^^ has been chosen with a monopole g^-dependence, np = 1, 
consistent with data in the D sector and the condition F^^ constant determines the parameter 
p^ — and then p^ — as a function of A^t^. A similar relation determines — and then A'^ — 
as a function of A2. 



25 



Ai{GeV) 


PL 


PL + ApL 




-Rj/* — ARj/-^ 


n2 = 2 


Ai, MAX = 8.112 


0.665 ±0.035 


0.700 


2.089 ±0.508 


1.581 


Ai = 6.810 


0.675 ± 0.033 


0.708 


2.263 ± 0.422 


1.841 


Ai, MIN = 5.426 


0.694 ±0.028 


0.722 


2.694 ± 0.494 


2.200 


n2 = I 


Ai, MAX = 6.113 


0.663 ±0.037 


0.700 


2.324 ±0.524 


1.800 


Ai = 5.770 


0.671 ±0.035 


0.705 


2.482 ± 0.475 


2.007 


Ai, MIN = 5.426 


0.681 ±0.033 


0.714 


2.714 ±0.514 


2.200 


n2 = 


Ai, MAX = 5.292 


0.660 ± 0.040 


0.700 


2.625 ±0.542 


2.083 


Ai = 5.237 


0.663 ±0.039 


0.702 


2.680 ±0.529 


2.151 


Ai, MIN = 5.183 


0.665 ±0.038 


0.703 


2.739 ±0.539 


2.200 



Table 3. 



The allowed domains for Ai,A2,Ay and Kp have been represented in 3 dimensional plots 
on Figures 12 to 17. However it might be useful to produce some numerical values obtained for 
Pl and Rj/qj in these domains and for that purpose, we have choosen, as an illustration, A2 = 
6 GeV , Ky = Kp = 'i) GeV, and for Ai, three values, Ai^ max, ^1, min and an intermediate 
value between the extremes. The results are represented in Table 3. A glance at Table 3 shows 
how difficult it is to fit simultaneously the large pp and the relatively small Rj/^, their opposite 
trends making the fit so difficult have been equally noticed [Q. 

The relative error on -Rj/^, is larger than the one of pi, and this feature is very useful for 
obtaining a fit. It is essentially due to the fact that Rj/isi, in addition to the errors on x^(0) 
and y^{0) that also enter in Api, has an uncertainty on z^{0) which is important. While the 
relative error on pj;^ is between 4% and 6%, the one of Rj/^ is between 18% and 24%. 

From these numerical results and those illustrated on Figures 12 to 17, it is clear that the 
scenario with a dipole form factor A2 is the one with the largest phase space in Ai, A2, A^, A^. 
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In this case it is relatively easy to accommodate both and Rj/^. We notice that the largest 
value of pl we can obtain in this model is pi = 0.694 ± 0.028 and for -Rj/ij, the smallest value 
is 2.089 ±0.508. 

For the scenario with a monopole form factor A2, the situation is less confortable even if it 
is possible to fit data on p^, and Rj/q,. Now the quantity pl varies by less than 3% and Rj/^ 
by about 17%. 

In the third case of a constant form factor A2, the allowed domain for the parameters 
Ai, A2, Av,Ap is very small and this possiblity, even if both p^ and Rj/q, satisfy the constraints, 
appears to be very marginal, the values of pi and Rj/q, being at the limits of the constraints. 

5°) The left-right asymmetry Alr for transverse polarization has not been experimentally 
measured. It has been defined in Eq. (pi.l8| ) and it depends only on the ratio y = y^{m?ji^) 
unambiguously related in our model to ^^(0). For our three selected models, we have 

The difference between the scenarios 77-2 = 2,1,0 is an allowed domain for Ai, A2 and Ay 
illustrated on Figures 12, 13 and 14. Scanning inside these domains, we make the predictions 

n2 = 2: 0.867 < Alr < 0.945 (V.15) 

722 = 1 : 0.837 < Alr < 0.910 (V.16) 

n2 = : 0.837 < Alr < 0.856 (V.17) 

The left-right asymmetry is large in the three selected cases, not as large as in the K* + ^i' 
case where it is close to one. We observe that the differences between the three scenarios are 
moderate. 



VI The decay modes B K[K*) + r}c 

1°) We start by considering the decay mode B — > K+rjc governed by the form factor Ff^^ [m^J. 
The Isgur-Wise relations determine F^^it^) from Fq^{{]) in a p^ dependent way : 

F^^{t°B) = Cbc [0.8460 - 0.0057 p^] F^^{0) (VI.l) 
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We notice that the coefficient of fi^ is very small as compared with the /i^-independent term 
in the bracket of Eq. (|VI.l| ). It follows that F^^{t°^) only weakly depends on fi^ . 

In Part II, we have choosen a monopole form factor for F^^{q^) and the parameter /i^ has 
been related to the pole mass Ai;' in such a way to obtain a form factor Fq^ independent on 
. From the previous considerations, the constant value of F^^ will be a weakly dependent 
function of Kp- With Kp in the range (5 - 6) GeV , F^^ increases slowly from 0.6286 ± 0.0251 



(for Ai. = 5 GeV) to 0.6327 ± 0.0253 (for Aj. = 6 GeV). The quantity Ff^(O) - . q 
been represented on Figure 6. 



BK 



(0) has 



2°) We now consider the decay mode B —>■ K* + rjc which is described by the hadronic form 



factors Aq^^* [m^^ 



^ The Isgur-Wise relation ( [1I.49| ) gives A^^* in terms of Af^ , A. 
the parameter A^. With the choice made in Part II for 



f^* and 



A^ 



m 



B 



m 



K* 



Ai 



(VI.2) 



A2 



A 



BK* , 



(VI.3) 



we get : 

^o^*imi) = 3-4602 Af^*{mlJ - 2.4602 

In Part V, we have obtained constraints on the scenarios allowed by the requirement to fit 
Pl and Rj/q,. They correspond to rii = —1 and n2 = 2, 1, 0. 

We introduce the function 5* (A) in order to relate the values of the hadronic form factors 
at = m^^ to the values at = t*£ where they are known from part II : 



S{A) 



A' 



'^B 



A2 



ni: 



and we get : 



A 



BK* 




3.4602 



[m. 



ric> 



S{A 



Af^*(t^°) - 2.4602 



Vc 



mz 



Al 



(VI.4) 



(VI.5) 



The form factor value Aq^ K'^ric' 



"2 Ms scenario- dependent, firstly by the value of ^2 ; ^2 = 
2, 1, 0, secondly by the values of the pole masses Ai and A2 in the restricted domains described 



in Part V. We vary Ai and A2 inside these domains and obtain : 



^2 
n2 



2 : 
1 : 
: 



A 



BK* 



m. 



A 



BK* 



A 



BK* 



m. 



[ni. 



1.3768 ±0.2928 
1.3979 ±0.2677 
1.3611 ±0.1790 



(VI.6) 
(VI.7) 
(VI.8) 
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We observe that the results for Aq^ (^^c) only weakly scenario- dependent due to cancella- 
tions between the two terms of Eg .( [VLSI) 

3°) The ratio i?^^ of these two decay modes 



is given by : 



R^^ = 0.3732 



(VI. 10) 



Using the values dVL^ ), ( |VrfD , ( [VLSD for ^^^'(m^J, and Eq.([VnD for F^^{ml), we obtain 

712 = 2 : Rvc = 1-7903 ± 0.7748 (VI.ll) 

n2 = l: Rrj, = 1.8456 ± 0.7221 (VI. 12) 

ria = : i?^, = 1.7497 ± 0.4810 (VI.13) 

Taking these results all together, we obtain the one standard deviation bounds 

1.02 < R^^ < 2.56 (VI.14) 

We comment, in the next Part VII, on the difference between these results and the bounds on 



Rr^^ previously obtained in Ref. ||T0| 



VII Comparison of Different Charmonium States Pro- 
duction 

1°) Ratio of decay widths with the same strange meson, K or K*, and different charmonium 
states are interesting quantities involving the leptonic decay constant fee- We define four such 
ratios referred to the most accurately measured J/"^ production : 

T{B ^ K + J/^) ' T{B ^ K* + J/<i/) ^ ^ ' 

^ ^ T(B^K + r^,) ^ YiB^K* + r^,) 

T{B ^ K + J/^) ' T{B ^ K* + J/^) ^ ■ ' 
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Assuming factorization and using the phase space estimates given in Part III, we obtain 



S* = 1.0057 



S = 0.4438 



T = 2.5180 



fr]c 



a - bxf + 2{l + c y ) 
(a - + 2(1 + c2?/2) 



T* = 0.8706 



2 



2 ^ 



(a - + 2(1 + cV) 



(VII.3) 

(VII.4) 
(VII.5) 
(VII.6) 



The rjc modes have not been experimentally observed and the only available information refers 
to J/\E' and modes. Using the PDG data ||5| collected in Table 1 of Part III, we obtain for 
S and S* values given in Table 4 : 



Ratio 


B+ 


B" 


B^, B° combined 


S 


0.68 ±0.32 


< 1.07 ±0.30 


0.68 ±0.32 


S* 


< 1.76 ±0.52 


0.89 ±0.59 


0.89 ±0.59 



Table 4 



2°) Using the values [y] of /j/* = (384±14) MeV and f^> = (282 ±14) MeV as estimated 
from the decays J/\E' e~^e~ and e^^e", we obtain : 

2 

I = 0.539 ± 0.066 (VII.7) 



and the quantity S is written : 

S = [0.2392 ± 0.0292] 



^,^^(m2 



J 1^1 



(VII.8) 



In our model, the hadronic form factor Ff^[q'^) has a monopole dependence with a pole 
mass Ap and we simply have : 

2 

(VII.9) 





2 




2 

- ^j/m 








-ml, 
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This ratio of form factors is a decreasing function of and so is the ratio S. At Aj? = 5 GeV, 
the prediction for 5* is : 

S{Af = 5 GeV) = 0.4363 ± 0.0537 (VII. 10) 

This prediction is in agreement, within one standard deviation, with the experimental value 
estimated in Table 2, Sexp = 0.68 ±0.32. Such an agreement continues to occur for larger values 
of A^ up to 6.27 GeV. 

The range of A^? depends on the three scenarios corresponding to n2 = 2, 1, and they are 
deduced from Figures 15, 16, 17 respectively. We get : 

n2 = 2: 0.4363 ±0.0537 > 5 > 0.3505 ± 0.0432 (VII.ll) 

n2 = l: 0.4363 ±0.0537 > S > 0.3790 ± 0.0467 (VII.12) 

n2 = 0: 0.4363 ±0.0537 > S > 0.4181 ± 0.0515 (VII.13) 

The errors quoted in Eq. (|VII.ll|) , ([VII. 12|) and ( |VII.13|) are due to the uncertainty on the 



leptonic decay constants f^' and fj/^. In conclusion, the theoretical predictions of our model 
for the three scenarios agree with experimental results within one standard deviation. 

3°) The analysis of the second ratio S* is more complex because of a large number of form 
factors involved. Using Eq.( |V11.7D , we get : 



S* = [0.5424 ± 0.0664] 



(a - hxf + 2(1 ± c2?/2) 



(VII.14) 



In our model the form factor y4f^*(g^) is linearly decreasing with a slope Ai, and we simply 
have 

Af^'(m^,) _ A^ - ml, 

We have computed the ratio 5** for the three scenarios n2 = 2, 1, using the values of Ai, A2 
and Ay inside the allowed domains obtained in Part V and represented respectively on Figures 
12, 13 and 14. 

The results of this scanning are : 

n2 = 2: 0.3287 ±0.0028 < S* < 0.4135 ± 0.0038 (VII.16) 

n2 = l: 0.3489 ±0.0034 < S* < 0.4015 ± 0.0039 (VII.17) 

n2 = 0: 0.3763 ±0.0039 < S* < 0.3867 ± 0.0040 (VII.18) 



31 



The errors quoted in Eqs.( |V11.16| ), (|V11.17D and ( |V11.18[ ) are computed in quadrature from 
those on the ratios f^' / fj/q,, x^(0) and y^{0). The theoretical predictions of our model for the 
three scenarios agree, within one standard deviation, with the experimental results estimated 
in Table 2 : S*^j, = 0.89 ± 0.59. 

4°) Kamal and Santra P] have studied the ratios S and S* denoted by them respectively 
as 1/R and 1/R . In the case of R, both monopole and dipole dependences for 

pBK 

are 

considered with a pole mass Ap = 5.43 GeV. Their conclusion is that a dipole behaviour for 
is needed in order to obtain an agreement for R between theory and experiment in the 
one standard deviation limit. 

The apparent contradiction between our result ( monopole for Ff^ ) and the one of Ref . 
is essentially due to the large experimental error of 47% for the quantity S or R. With 6 = 0.47 
the relation at first order in 6, (1 ± = 1 =F 5 is not valid and one standard deviation limit 
for S and one standard deviation limit for R are different concepts. However, since the main 
part of the experimental error is due to the K + '^' mode and for that reason the consideration 
of one standard deviation for S ( where K + enters in the numerator ) seems to be more 
relevant than for R. 

A similar situation occurs for S* and R' . Here the experimental error is even larger, 66.7%, 
and it is mainly due to the K* + \E'' mode which enters in the numerator of S*. Again the 
one standard deviation limit for S* and the one standard deviation limit for R' are different 
quantities. 

Also the pole masses in Ref.0 are taken only at some fixed values, while in our approach 
these poles sweep inside the allowed domains of Figs. 12 — 17. 

For the ratio R' as previously done for p^, they propose seven scenarios. Furthermore, con- 
sidering only in the one standard deviation limit for R' , they exclude four scenarios where Af^* 
is either constant or linearly decreasing with and conclude that if factorization assumption 
were to be held, then the only scenarios that are consistent with experiment are those in which 
rises with g^. We observe however that R' (or S*) is not an independent ratio but re- 
lated to the other ratios by S*Rj/q, = SR^', such that considering R' (or S*) alone might be 
inadequate. 

5°) Comparing the K + rj^ and K + J /'^ decay modes, we now consider the ratio T which de- 
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pends on the ratio of the decay constants /^^ and /j/*. Unfortunately /^^ is not experimentally 
known and we use theoretical estimates if we want to make predictions. 
Rewriting Eq (|VII.5|) in the form : 

T = 2.5180 ( 4^] SviAp) 



where 



m 



m 



(VII.19) 



(VII.20) 



We compute Sy in our model where F^^^ is constant and Ff^ has a monopole dependence 
with the pole mass Ap. As a consequence, we simply have 

.2 \ 2 



1 



m 



J/* 



A2 



(VII.21) 



The function Sy is an increasing function of A^?. The allowed values of Ap have been discussed 
in Part V (Figs. 15 — 17), and we obtain : 



ria = 2 : 5 GeV < Ap < 5.71 GeV 
n2 = l: 5 GeV < Ap < 5.39 GeV 
n2 = : 5 GeV < Ap < 5.10 GeV 



, 0.38 < Sv < 0.50 (VII.22) 
, 0.38 < Sv < 0.45 (VII.23) 
, 0.38 < Sv < 0.40 (VII.24) 

As pointed out in Ref . Ijl^ , a measurement of the ratio T will provide an opportunity to extract 
the scalar decay constant /^^ from experiment. Unfortunately the decay mode B ^ K + 7]^ has 
not been experimentaly observed. 

There exist various theoretical ways to estimate both f^^ and /j/<t. Using the estimate 



quoted in Ref. 1 17 



frjc 



0.993 



we make predictions for the ratio T : 

n2 = 2: 
n2 = l : 
n2 = 0: 



0.94 < T < 1.24 
0.94 < T < 1.12 
0.94 < T < 0.99 



(VII.25) 

(VII.26) 
(VII.27) 
(VII.28) 



33 



6°) We finally discuss the last ratio T* which, under the factorization assumption, has the 
form given in Eq.( [VII.6D . We compute T* for the three scenarios n2 = 2, 1, 0, using the values 
of Ai, A2 and Ay inside the allowed domains obtained in Part V. By writing 



we obtain 



Vc 



r 



(VII.29) 



n2 = 2 : 
^2 = 1 : 
^2 = : 



0.6148 ±0.1108 < T < 0.7740 ±0.1002 
0.6097 ±0.1259 < T* < 0.7375 ± 0.1347 
0.6015 ±0.1415 < T* < 0.6217 ±0.1444 



(VII.30) 
(VII.31) 
(VII.32) 



If one accept the value ([VII. 2 5|) for the ratio of leptonic constants, we make predictions for the 
ratio T* which are weakly scenario-dependent : 

(VII.33) 
(VII.34) 
(VII.35) 



n2 


= 2 : 


0.50 


< 




< 


0.86 


n2 


= 1 : 


0.48 


< 




< 


0.86 


n2 


= : 


0.45 


< 




< 


0.76 



7°) These results are now compared with those obtained in a previous paper ||T0[. The 
scenario- dependent parameter Sy defined by Eq. (|VII.20D is obviously very different from its 
value in the BSW model considered in Ref. [|1(J|]. The second scenario dependent parameter Sa 
entering in T* : 



A 



BK* , 



2 ^ 



aBK* 
^0 



(0) 



Af^*(0) 



(VII.36) 



is computed in our model with the result : 



n2 = 2 : 

^2 = 1 : 
7^2 = : 



0.767 < Sa < 0.793 
0.832 < Sa < 0.849 
0.966 < Sa < 0.978 



(VII.37) 
(VII.38) 
(VII.39) 



For the BSW model H, is 1.069 . 
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The ratio Sa/ Sy entering in the ratio i?^^ is more dependent on the scenario than separately 



in 


our model 


are : 












= 2 : 


1.53 


< 


Sa/ Sy 


< 


2.09 




= 1 : 


1.85 


< 


Sa/Sv 


< 


2.19 


n2 


= : 


2.44 


< 


Sa/Sv 


< 


2.54 



(V1I.40) 
(V11.41) 
(V11.42) 



The bounds used in our previous paper [|l^, 1 < Sa/Sv < 1.4 are largely underestimated, 
essentially because of the behaviour in of the hadronic form factors : F^^ constant and 
Af^ linearly decreasing with . Therefore, the predictions of Ref. |]T0| for i?^^ are smaller 
than those obtained here. 



VIII Radiative Decay B ^ K* + 

1°) The radiative decay B K* + 7 does not occur at the tree level in the standard model. 
At the one loop level, we have the so-called Penguin diagrams and, for the case consider here, 
the dominant contribution is the one corresponding to the exchange of a virtual t quark. In 
this approximation, the branching ratio is given by : 

BRiB^K* + j)= f ^ \V4' \V,s\' ml (l-^X \Cr{m,)\^ { } (VIlI.l) 

h 256 TT'^ \ ITLg J 

where the quantity { } depends on the hadronic form factors associated to the tensor and 
pseudotensor current taken at = and for real photon [l^ : 



{ }=(l + -)' \v^m'+{l-—\ I4(0)P (VIII.2) 
\ nib/ \ nib/ 

The quantity C-j{nib) is the Wilson coefficient associated to the relevant weak current. It takes 

into account the large QCD corrections |^ and it plays a determinant role in the numerical 

calculation of the rate. 

We use the unitarity relation of the CKM matrix 

Vtb = -v,b - v:, (VIII.3) 



and we neglect the second term in the right hand side of Eq. ( |VIIL3D being 0{sin'^6c) with 
respect to the other ones. 
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Numerically we take 

\V^s\ = 0.970 , \Va\ = 0.040 , tb = 1.54 ■ 10"^^ s , = 4.7 GeV (VIII.4) 



and for the Wilson coefficient, choosing rrit = 174 GeV and Aqcd = 200 MeV, we have |]19 
Cjirrih) = 0.325. 
The result is 

BR{B ^K* + -f)= 4.45 ■ 10"^ { } (V111.5) 



The radiative decay mode B K* + 7 has been experimentally observed by CLEO ||20| 
BR{B ir* + 7) = (4.5 ± 1.5 ± 0.9) ■ 10"^ = (4.50 ± 1.75) ■ 10"^ (VIII.6) 
and we get 

{ U = 1-011 ±0.393 (VIII.7) 

2°) Following Isgur and Wise we assume that in the B meson at rest, the b quark 
spinor in the weak current satisfies the free Dirac equation for a spinor at rest : 706 = b. As a 
consequence, the tensor (pseudotensor) current is related to the vector (pseudovector) current 

Q bo, lj]b = - 2q-fjb , q [7^, 7^-] 75 6 = 2 g 7^- 75 6 (VIII.8) 

It is then straight forward to compute the four tensor or pseudotensor hadronic form factors 
in terms of the ususal BSW form factors V,Ai,A2 and ^3. The two relations of interest here 
are : 

t;^(0) = a^(0) = f 1 + Af^*(0)+fl-^') V^^\0) (V111.9) 

In the limit f"^(0) = 02 (0), the quantity { } of Eg. ( VIII. 2 ) becomes in the approximation 



{ } = 2K;^(0)p (VIII.IO) 
and using the estimate ([VIII. 7| ) for { }, we obtain 

^"^(0)^.^ = 0.71 ±0.14 (VIII.ll) 

The tensor and pseudotensor hadronic form factors have been computed in various models, 
quark constituents models ||2T|, vector meson dominance models [^, lattice gauge theories 
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||23|| . In most of the estimates, the equahty f"^(0) = a^(0) is obtained with v'^{0) in the range 
0.5 to 1, e.g., consistent with the CLEO resuh ( |VIII.11|) . A recent estimate of Griffin, Masip 
and Mc Guigan using the Isgur-Wise relation gives f^(0) = 0.97 ± 0.13 . 



It is amusing to observe that the BSW model ||^ which fails [|I| in explaining the ratio pi 
and Rj/^ in B K{K*) + J/^ produces for f^(0), using Eq.( |VIII.9|) , a value f^(0) = 0.69 in 
very nice agreement with experiment. It is dear that the radiative decay B — > K* + 7 is not a 
very efficient filter for models. 

3°) In our model, we compute f^(0) = aj(0) using formule ( VIII. 9 ). For Af^ we use a 



linearly decreasing function of and for V^^* a increasing function of of the dipole type, 
the input values being predicted at = t*^ by the Isgur-Wise relations. As a consequence, we 
obtain a relatively large Af^*{0) and a relatively small V^^*{0). The minimal value of v^{0) 
corresponds, in the allowed domain of the parameter space, to A2 = 6 GeV, Ay = 5 GeV, 
Ai = Ai^ MAxi-^v, -^2) and the results are in the one standard deviation limit, including QCD 
corrections : 

(VIII.12) 
(VIII. 13) 
(VIII. 14) 

If the assumption v'^{0) = a.^(0) is correct, we see that the scenario ^2 = 2 can accomodate 
the experimental result ([VIII. 11|) . A fit is clearly more difficult for the scenario n2 = 1 and it 
seems to be impossible for the scenario n2 = 0. 

However we must be aware of the fact that the estimate of the QCD correction, which is 
scale dependent, may have some uncertainty which has been disregarded in the experimental 
error quoted in Eq. ( |VIII.11| ). A theoretical error has to be added which might be as large as 

15 % m. 



^2 


= 2 : 


.^(0) 


> 


0.94 




= 1 : 




> 


1.11 


n2 


= : 


v^{0) 


> 


1.32 



IX D K(K*) Hadronic Form Factors 



1°) The B —>■ K{K*) and D K{K*) hadronic form factors are related by the SU(2) heavy 
flavour symmetry of Isgur-Wise. From the considerations of Part II, it is clear that the g^- 
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dependence for the form factors Fi, Ai, A2 and V, namely 

(i) Same values for rii, 77-2, ny and np in both B and D sectors; 

(ii) The pole masses in these sectors are related by Eqs.( II.2"^) and (P^I.37|) of Part II. 
For Fq and Aq the situation might be different in the two sectors but the g^-dependence of 

Fo is known from that of Fi and for Ao, its g^-dependence follows from that of Ai and A2 as 
explained in Part II. 

2°) Let us first consider the semi-leptonic decay of D mesons. The relevant hadronic form 
factors are Ff^ hi D ^K + i+ + ue and Af^\ A^^\ V^^' for D ^K* + i+ + u^. Using 
the dimensionless variable t = q'^/mj), we introduce the dimensionless function X{t) : 

A-(*) ^ i § (.X.1) 

which is independent of all parameters entering in the semi-leptonic relevant scale factor. 

We recall that the quantities x^(0), y^{0) and 2:^(0) (used in this paper for normalizing 
the B sector) have been extracted from experimental data on semi-leptonic decay in a scenario- 
dependent way, because the variation with of the form factors F^^ , Af^*, A^^* and V^^* 
has not been measured. 

3°) The distribution in semi-leptonic decay D ^ K + + ui is written in terms of the 
hadronic form factor F^^{q^) as 

+ = ^ \VMK{q')r \Fri<l')\' (IX.2) 

where the dependent momentum K{q^) is given by 

K{q') = ^ {{ml + ml- q'f - (IX.3) 

In the zero lepton mass limit, < g^ < [vtid — vtlk)^ . 

In our model Ff^{q'^) has a monopole g^-dependence with a pole mass Ajjp related to 
in the B sector by Eq. ( |II.24|) . Defining the dimensionless parameters : 



^DF 



aF = jf- (IX.4) 



we obtain for X{t) the expression : 



I [aF) (1 — ap ty 
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where the integral /(a^) is defined by the normahzation condition X{t) 



{1 — ap x) 



I{ap) = I / " ^ dx (IX.6) 



Of course, the semi-leptonic rate is simply given by : 

TiD^K + i^ + u,) = ^^ \V^f |i^r(0)P /(M (IX.7) 

The normalized distribution X{t) for the semi-leptonic decay mode D K + + ug is 
represented on Figure 18 for values of ap corresponding to the bounds on A^? obtained in Part 
V and illustrated on Figures 15, 16 and 17. 

na = 2 : 5 GeV < Ap < 5.71 GeV 1.097 > > 0.630 (IX.8) 
n2 = 1 : 5 GeV < Ap < 5.39 GeV 1.097 > ap > 0.787 (IX.9) 
n2 = 0: 5 GeV < Ap < 5.10 GeV" 1.097 > > 0.999 (IX.IO) 

The distribution X{t) is a monotonically decreasing function of t. Its shape is not very sensitive 
to ap except in the neighbourhood of t = 0. 

An estimate for the slope of the distribution at = has been given by Witherell W2 



using two models for the g^-dependence of Ff^lq"^), an exponential form and a monopole form. 
The result consistent for the two models translated in the ap language is 

0.76 < ap < 1.30 (IX.ll) 

which means in our model 

5.36 GeV > Ap > 5.02 GeV (IX. 12) 
Therefore our bounds ( [IX. 8| ) — ( |IX.10| ) are consistent with experiment. 



4°) The distribution in the semi-leptonic decay D — > K* + + depends on the three 
hadronic form factors A^^ (g^), (g^) and V^^* [q^). We have three possible polarizations 
for the final K\ X = 0,±1. 

It is convenient to define dimensionless parameters : 

r* = ^ , t = \ , a, = -f- ■ j = l,2,V IX.13 

rriD mf) Ajjj 
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where the pole masses in the B and D sectors are related by Eg. ( [11.3'^ ) . 
The fixed distribution is given by : 



dT{D + ue)x Glm 



2 ^5 
D 



dt 



1927r3 



il + r*y \A^^ MA(t;ai,a2,ay) (IX.14) 



where 

Mi(t;ai,a2) = k{t) 

M±{t]ai,av) = t k{t) 
with k(t) given by : 



a{t) 



(1 - aity 
1 



Tcit) 



(1 - a2t) 



(1 - (1 - avty 

1/2 



k{t) = {(l + r*2_t)2_4r*2| 
The coefficients d{t), b{t) and c{t) are : 



a(t) 

m 
m 



-t 



2r* 



2r*(l + r*y 
kit) 



[l + r 



*\2 



We define the integrals : 



Ix(ai,a2,av) 



and the integrated rate is given by : 



(l-r-)2 



Mx{t;ai,a2,av) dt 



1927r3 

The normalized distribution is : 

X{t) 



IK:, 



Sa Mx{t;ai,a2,av) 



(IX. 15) 
(IX.16) 

(IX. 17) 

(IX.18) 
(IX.19) 

(IX.20) 



(IX.21) 



Sa Ix{ai,a2,av) (IX.22) 



(IX.23) 



Sa Ix{t;ai,a2,av) 

Of course for the scenarios selected in Part V, we have rii = ~l,ny = 2 and three possible 
values for n2', ^2 = 2, 1,0. We have computed X{t) in these three cases by using the PDG 
values for x'^{0) and y^{0). The parameters ai,a2,Q:v — or equivalently Ai,A2,Ay — are 
constrainted to stay inside the allowed domains represented on Figures 12, 13 and 14. The 
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results are shown on Figures 19, 20 and 21. As in the previous case, the largest sensitivity of 
X{t) to the parameters is in the neighbourhood of t = 0. 

In an analogous way, it is possible to study the distributions for the polarization param- 
eters pf{t) and Afji{t) : 



La Mx{t;ai,a2,av) 
^si ^ M_(t; ai, q;2, ay) - M+{t; ai, ^2, ay) (1X25) 
M_{t;ai,a2,ay) + M+{t;ai,a2,ay) 

We only give here the integrated ratios pf; and Af^j^ where the functions Mx{t;ai,a2,ay) in 
Eg. ( [IX. 151 ) and ( [1X.16| ) are replaced by their integrals over t, Ix{ai, a2, ay). The results for the 



three cases r;,2 = 2, 1, are the following 



712 = 2 : 0.516 < pi < 0.541 (IX.26) 



0.885 > Afji > 0.829 



n2 = 1 : 0.526 < pf < 0.541 (IX.27) 



0.904 > Aiji > 0.857 



n2 = : 0.536 < pf < 0.538 (IX.28) 

0.904 > Ai^ > 0.892 

In Eqs. (|IX.26| ) - ([1X.28|) the results are presented in such a way to exhibit a correlation between 
the largest (smallest) pf and the smallest (largest) Afj^. 

5°) Let us now consider the hadronic two body decays D° K{K^*) + 7r+(p+). Assuming 
factorization to be justified in the D sector, the various decay rates can be determined by using 
the hadronic D K(K*) form factors studied here. 

At the quark level, the tree level diagram is of the spectator type and the various rates are 
written in the following way : 

BR = BRo ■ ( ^\ ■ PS ■ FF (IX.29) 

\mnj 



41 



The scale BRq is given by : 

BRq 



I V2 ) 



\ycs\ IVudl — — 0-1 



Sir 



h 



(IX.30) 



where Oi is the BSW [|] parameter for color favoured processes. As in the Part III, PS is a 
phase space factor, /^^ is a leptonic decay constant /^+ or /p+ which are experimentally known 
and FF depends on the hadronic form factors. 
The results are : 



a) D" - 


R- 


+ vr+ , 


PS = 


0.3993 , 


FF = 




{ml)? 


b) D" - 


R- 


+ , 


PS = 


0.1936 , 


FF = 


iFn 


ml)? 


c) D° - 


^ R* 


' + 7r+ 


PS = 


0.2216 , 


FF = 




{ml)? 


d) D" - 


R* 




PS = 


0.0843 









(IX.31) 
(IX.32) 
(IX.33) 



FF = |Af^*(mJ)p {{an - boXo)' + 2[1 + coW]} (IX.34) 



where 



and 



an = 1.5310 



bo = 0.2419 



A?^*{ml] 



CD = 0.2087 



Vd 



(IX.35) 



(IX.36) 



In order to test the assumption of factorization in the D sector, we shall study two ratios of 



rates : 



T{D° R*~ + 7r^ 
T{D" ^ R- + n+ 



T, 



Using Eq.( pC3lD and ( pO^ ), we obtain 

R = 0.5550 



DK* 



D 



.ml 



T{D° R- +P+) 
r(F)° ^ ir- +7r+) 



(IX.37) 



mt 



(IX.38) 



It is legitimate to neglect the variation of the form factors Fq and Aq between = and 
= ml. Using the normalization conditions M, 



Fo^^(O) 



aDK' 
^0 



(0) 



F/^^(0) 
mD + rriK* 



2m 



K* 



(0) 



mo — mx* 



2m 



K* 



A^' 



(0) 



(IX.39) 
(IX.40) 
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and the values at = of the form factors as given by the PDG : 

F/^^(0) = 0.75 ±0.03 , Af^*(0) = 0.56 ±0.04 , ^*(0) = 0.40 ± 0.08 (1X.41) 
we obtain : 



A^^*(0) = 0.6473 ±0.0757 
and the theoretical prediction for i?^ is 

(i?^)^^ = 0.4134 ±0.1022 

The experimental result is : 



(IX.42) 



(IX.43) 



(IX.44) 



(i?^)e..p = 1.22 ±0.16 

The discrepancy between theory and experiment is very large. The theoretical prediction is 
clearly scenario-independent and it is unlikely that final state interaction would be able to 
fill the gap between theory and experiment. The most probable explanation is a failure of 
factorization in the D sector. 

Let us now consider the quantity T^. Using Eg. ( |IX.31|) and ( PX.32|) , we get 

2 

Tz) = 0.4850 (^1 T;W7rf^ (IX.45) 



For the leptonic decay constants /^r and fp, we use the experimental values 

= 131.7 MeV fp = 212 MeV 

and Td is written as : 



Td = 1.2566 



(IX.46) 



(IX.47) 



Neglecting as previously the variation of Fq between = and = ml and using the 
normalization condition (p!X.39|) , we get 



F^^{ml) F^^{ml) 



F^^{ml) - Fi^^(O) 
In this paper we have used a monopole g^-dependence for Fi, n^? = 1 and we have 

^DK 



(IX.48) 



i^r(o) 1 



(IX.49) 
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where the pole masses and A£,f in the B and D sectors are related by Eg. ( [11.241) . With the 
constraint > 5 GeV, the maximal value of is obtained with A^p = 1.78 GeV — which 
corresponds to Ai;' = 5 GeV — with the result. 

(Td)^^ < 1.8946 (IX.50) 

The experimental value is significantly larger : 

(Td),,, = 2.59 ± 0.34 (IX.51) 

Again the factorization assumption seems to be in a bad shape in the D sector. 

X Critical Discussions and Conclusions 

1°) We have shown that the assumption of factorization is not ruled out by experimental data for 
the colour-suppressed decay modes of the B meson, B + J/^(^ ). The failure pointed 

out in Ref. |l[] might be due to inadequate choices of hadronic form factors and the aim of this 
paper is essentially to exhibit possible g^-dependences that are able to explain experimental data 
and particularly the ratios pi and Rj/qj. Of course the possibility to understand experiment is 
not necessary a proof of factorization. 

Let us first summarize the assumptions and constraints contained in our model. 

(A) Assumptions : 

1. Factorization holds for color supressed B decays and final state strong interaction effects 
can be neglected. 

2. The SU(2) heavy fiavour symmetry between the b and c quarks is realized by the Isgur- 
Wise relations 0. 

3. The input experimental data in the D sector are taken from the analysis of semi-leptonic 
decays D ^ K + i'^ + U£ and D —>■ K* + + in the form of the values at = of the 
D —>■ K{K*) hadronic form factors^]. 

(B) The experimental constraints are : 
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1. The experimental rates for B K + J/"^, B K* + J/'^, B K + ^' and B ^ K* + '^' 
used in the form of ratios of rates Rj/^, R^' , S and S* [defined respectively by Eq. ([III.19| ) 
and Eq. (|VinD ]. 



2. The observed longitudinal polarization fraction piin B ^ K*+J/'^ [defined by Eq.( |III.l"5D ] 
(C) The theoretical constraints are : 

1. The explicit form of the dependence of the hadronic form factors Fi, Ai, y42, V choosen 
as [1 - gVA^]~" with n = -1, 0, 1, 2. 

2. The pole masses A of the various form factors in the B sector are limited to the (5—6) GeV 
range in order to relate them in a likely way to bs bound state masses. 

3. The ratios of form factors ^i^{q^) and A^(g^) defined in Eqs. ( [11. 1^ and ( |11.46| ) are assumed 
to be independent of and related in a natural way to the pole masses Kp and A2 by 
Eqs.dg) and (|ir55| ). 

2°) Among our three theoretical assumptions, the weakest one seems to us the third, i.e., the 
use of the experimental values at = of the D K{K*) hadronic form factors as deduced 
in a scenario-dependent way from the experimental triple angular distribution of semi-leptonic 
decay. We must emphasize that measurements of the g^-dependence of these form factors are 
not available and their values at = are obtained by extrapolation at g^ = of experimental 
data, assuming monopole g^-dependence for all form factors. 

Consider first the simplest case of the form factor FP^{q^). The two decay modes 
K i+ ue and D" K~ i+ Ui are expected to have equal rates, the associated weak current 
cs being isoscalar. Experimentally these two rates differ by few standard deviations [^ and 
some average dominated by the most accurately measured D° mode is used for extracting the 
quantity F[^^{0) and the slope of F['^{q^) at g^ = from the measured g^ distribution. In 
order to have an estimate of the uncertainties of the analysis, let us notice that the rate value 
used by Witherell [0] is 10 % higher than the one of the PDG |^ both values being given with 



5 % errors. In spite of this difference, the values obtained by the PDG [§, F/^^(0) = 0.75±0.03 
and by Witherell [0, F/^^(0) = 0.77 ± 0.04 are very similar. 
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We have considered the problem of determining Ff^{0) from the semi-leptonic rate by 
using Eq. ([[X.7| ). The resuh depends on the parameter ap simply related to the pole masses 



Adf and Ap in the D and B sectors. 
The results are : 

n2 = 2 : 5 GeV < Ap < 5.71 GeV , 0.68 < F^'^{0) < 0.79 (X.l) 

712 = 1: 5 GeV < Ap < 5.39 GeV , 0.68 < ^(0) < 0.77 (X.2) 

n2 = : 5 GeV < Ap < 5.10 GeV , 0.68 < Ff'^iO) < 0.73 (X.3) 

The numerical value of the rate used in this analysis is the PDG value of (8.2 ± 0.4) 10^'^ 
We notice that the slope results presented by Witherell ||12| correspond in our language to 
5.02 GeV <Ap< 5.36 GeV. 

We now comment on D ^ K* Vf data, for which the g^-dependences for the relevant 
hadronic form factors, A^* , A^^' and V^^* , are not yet available. By fitting the distribution 
and the angular distribution due to the K* ^ K + n decay, the ratios x^i^Q) and y^{Q) 
are obtained assuming for the form factors a monopole g^-dependence with the pole masses 
given by the lowest lying cs bound states, -D*(2110) and D**(2560). The results of the PDG 
and of Witherell are nearly identical. The quantity Af^*(0) is then obtained from 



the transition rate. In this case the semi-leptonic decay rates T{D° K*^ + + vi) and 
r(D+ — > K*" + £^ + Vi) are equal within errors. However the average value of the rate used by 
the PDG and by Witherell differs by 11% and as a consequence the values of y4f^*(0) obtained 
in the two cases are different : 



PDG[h\ Af^'(O) = 0.56 ±0.04 , ^'(0) = 0.61 ± 0.05 (X.4) 



We study the problem of determining A^^ (0) from the semi-leptonic decay rate by using 
Eq.( [lX.22D . We vary the parameters Ai, A2, Ay (ai, 02, ay) in the allowed domains determined 
in Part V and the results are : 

n2 = 2: 0.601 < Af-^*(0) < 0.645 (X.5) 

n2 = l: 0.618 < Af^*(0) < 0.645 (X.6) 

^2 = 0: 0.638 < Af^*(0) < 0.647 (X.7) 
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Our Af^ (0) is to be compared to the PDG one 0.56 ± 0.04, since the numerical value of the 
rate used in our analysis is the PDG value (4.6 ± 0.4) 10^° s^^. 

The errors quoted for Ff^{0) and Af^ (0) in references and [12| do not take into account 
the theoretical uncertainties due to the g^-dependence used for analysing the experimental data. 
The formal results previously obtained with our model illustrate clearly the point and it seems 
to us that the errors have been underestimated. 

Finally let us point out that the various ratios studied in this paper have different types of 
dependence with respect to the form factor values at = in the D sector. 

(i) . Alr and depend only on ^^(0). 

(ii) . pl and depend on x^(0) and y^{0). 

(iii) . Rj/^ and Rq,' depend on x^(0), y^{0) and z^{0). 

(iv) . i?^^ depends on x^(0) and z^{0). 

(v). S and T are independent of these three ratios, 
(vi). S* and T* depend on x^(0) and y^{0). 
For the semi-leptonic normalized distribution X{t), it is independent on these ratios in the 
D ^ K + + case and it depends on a;'^(0) and y^{0) in the D K* + + i/£ mode. 

3°) In addition to the necessity of an improvement in accuracy for the observed rates — in 
particular, those involving the ^' are still badly known with errors at the 45% or 64% level — 
we need a detection of the K + 7]^ and K* + rjc modes. On the polarization side, only the ratio 
Pl has been measured and it was an important quantity to be taken into account for the fits. 
It is clear that the knowledge of Alr, Pl ^^"^ -^lr predicted by our model would considerably 
help in reducing the size of the allowed domains in the parameter space. 

The best situation would be to select only one scenario and a small domain in the parameter 
space. The worse situation for our model would be that the new measurements exclude the three 
presently remaining scenarios. Our approach being purely phenomenological is not connected 
with any theoretical calculation like, for instance, QCD sum rule or lattice gauge theories and 
our model is certainly not the unique way to compute B K{K*) hadronic form factors. 

However if we are in the best situation previously mentioned, it will be necessary to provide 
a theoretical support to the so determined hadronic form factors and for that, results of Ref.p| 
seem to be in a good shape because of the unusual behaviour prediction for Ai{q'^). 
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If we are in the worse situation, it will be reasonable to think seriously of the role played 
by non-factorizable contributions. 
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Appendix A 

1. Error on pi ( Eq: (|ITLT7D ). 



ApL = 2pL{l- Pl) 



hx 



a — bx 



2 / A D \ 2 



X 



D 



,2„.2 



2c'y 



1 + c^y'^ 



'Ay 



1/2 



(A.l) 



2. Error on = (a - bxf + 2(1 + c^y'^^ 



AA^ = 2 <^ [bx{a - bx)f 



'Ax 



D\ 2 



X 



D 



+ 



2cV 



2 / Ay 



D\ 2' 



1/2 



(A.2) 



3. Error on Rj/^ ( Eq: (|nLT9D ). 



ARjjm — Rj/m 



1/2 



(A.3) 



4. From the PDG values in the D sector at = 



Ax 



D 



X 



D 



15 
73 



Ay 



D 



25 
189 



'Az 



D\ 2 



Ti) + Q' 



(A.4) 



5. In the ^' case, we make the situations a ^ a\b ^ b' , c ^ c , x ^ x\ y y' and z ^ z . 



6. Error on S ( Eq:( [Vir3D ). 



AS = 2 S 



'Afj/A , (Af^ 



2\ 1/2 



+ 



(A.5) 
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7. Error on S* ( Eq:(|VITl]) 



2 S* 
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'2S 



1 2 



+ 



bx{a — bx) b'x'{a' — b'x' 



-\ 2 



N 



+ 4 



cV 
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AT' 



i2 i2 



Ax 



D\ 2 



X 



D 



c y 
N' 



'Ay 



1/2 



(A.6) 



Appendix B 

Semi-leptonic decay D ^ K Ve '■ 

Using for the hadronic form factor F^^{q^) a monopole g^-dependence, the decay rate is 
given by Eq.( |lJK-7| ) where the integral /(a^) is defined in Eg. ( |1X.6] ). From the experimental 
rate, it is possible to deduce the value of Ff^^{0) if the slope parameter ap is known. 

The computation of Ilap) has been done both numerically and analytically. In the case of 
interest here : 



the result is 



1 



< ap < 



1 



rriD 



(B.l) 



I(aF) 



+ 



1-r' 

3 

oi 



{-6 + 9(1 + r^)aF - 2(1 - r^fal] 



-J {l-2(l + r2)a^ + (l + r^)a^}lni 
1- 



F 



)aF\J[{l + rYaF-l\[l-{l-rfaF] 



- -Arctg ; ^ ^ ' = } (B.2) 

2 (1 - r2)y^[(l + ryaF - 1] [1 - (1 - r)^^^] J 

Of course such a formula can be used for any semi-leptonic decay of the type 0^ — > 0^ + 
+ Ue provided the corresponding form factor has a monopole g^-dependence. 
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Figure Captions and Table Captions 
A. Figure captions 

1. Figure 1 : Parameters for the B ^ K transitions. On the horizontal axis is A^., on the 
left vertical axis is the parameter /x, on the right vertical axis is the pole mass A|,p. The 
relevant relations are 

B _ ml - m\ _ fi^ -a _^ _ m]j-m\ 

Adf = — Ap- [1 [mbmc - rrij^) 

rub \ mbJ 

The values of A^ and Aj:,p used by Bauer-Stech-Wirbel are indicated by the point BSW 
and they are A^ = 5.43 GeV and Adf = 2.11 GeV . 

2. Figure 2 : Parameters for the B K* transitions. On the horizontal axis is Ag, on the 
left vertical axis is the parameter A, on the right vertical axis is the pole mass A|)2. The 
relevant relations are 

A — TT, , A — , A 



Ai ' 1 + aA^ ' Al^ 

Ad2 = — - ( 1 - — ) {mbrric - m^*) 
rUb \ rUbJ 



rUb \ mb. 

The values of A2 and A/)2 used by Bauer-Stech-Wirbel are indicated by the point BSW 
and they are A2 = 5.82 GeV and Ad2 = 2.53 GeV. 

3. Figure 3 : The quark level colour-suppressed diagram bq —>■ (cc) -|- {sq). 

4. Figure 4 : Af^*(m^,) and x^{ml,) as functions of A2 for 5 GeV < A2 < 6 GeV. One 
standard deviations are indicated with dotted points. 

5. Figure 5 : The fractional longitudinal polarization as a function of A2 for 5 GeV < 
A2 < 6 GeV. One standard deviations are indicated with dotted points. Our scenario- 
independent upper limit < 0.5664 is indicated. 
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6. Figure 6 : The hadronic form factor Ff^{q'^) at = 0, = mj^^ and = rn^, as 
functions of A^? for 5 GeV < Ap < 6 GeV. One standard deviations are indicated with 
dotted points. 

7. Figure 7 : The ratio of rates R^' as a function of A2 and Ap for 5 GeV < A2,Af < 
6 GeV. The cross points indicate the one standard deviation experimental hmits 0.44 < 
R^' < 3.62. 

8. Figure 8 : The allowed domain in the Ai, A2, Ay space due to the constraints p^ + Ap^ > 
0.70 for A2,Av G (5 - 6) GeV, Ai > 5 GeV. The scenario is [ni,n2,nv] = [-1,2,2]. 

9. Figure 9 : Same as Figure 8 for the scenario [-1, 1, 2]. 

10. Figure 10 : Same as Figure 8 for the scenario [-1, 0, 2]. 

11. Figure 11 : Same as Figure 8 for the scenario [-1, 2, 1]. 

12. Figure 12 : The allowed domain in the Ai, A2, Ay space due to the constraints pl + 
ApL > 0.70 and Rj/^ - ARj/^ < 2.20 for A2,Ay e (5 - 6) GeV, Ai > 5 GeV. The 
scenario is [711,^2, riy] = [—1,2,2]. Ap^ and ARjj^^ are theoretical errors induced by 
experimental errors of x^{Q),y^{Q) and ^^(0). 

13. Figure 13 : Same as Figure 12 for the scenario [-1, 1, 2]. 

14. Figure 14 : Same as Figure 12 for the scenario [-1, 0, 2]. 

15. Figure 15 : The allowed domain in the Ap, A2, Ay space due to the constraint 

Ai, min{A2,Av,Ap = 5 GeV) < Ai < Ai, max(A2, Ay) with A^ > 5 GeV. The scenario 
is [ni,n2,ny] = [-1,2,2]. 

16. Figure 16 : Same as Figure 15 for the scenario [-1, 1, 2]. 

17. Figure 17 : Same as Figure 15 for the scenario [-1, 0, 2]. 

18. Figure 18 : The normalized dimensionlcss distribution X{t) for the semi-leptonic decay 
D ^ K + l'^ + The scenario n2 = 2 corresponds to 5 GeV < Ap < 5.71 GeV , the 
scenarios 712 = 1 corresponds to 5 GeV < Ap < 5.39 GeV and the scenario ^2 = 
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corresponds to 5 GeV < Ap < 5.10 GeV. By Eg. ( [11.24]) the pole masses A^^;' in the 
D sector can be obtained from Ap given here. 

19. Figure 19 : The normahzed dimensionless distribution X{t) for the semi-leptonic decay 
D — > K* + + Ui. The scenario is [rii, -0.2, ny] = [—1, 2, 2]. The thickness of the curve is 
due to the Ai, A2, Ay ranges. 

20. Figure 20 : Same as Figure 19 for the scenario [-1, 1, 2]. 

21. Figure 21 : Same as Figure 19 for the scenario [-1, 0, 2]. 
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B. Table Captions 

1. Table 1 : 

Experimental data of the decay rates for K^*^ + \E'' and i^T*^*^ + J/^ as averaged by PDG 

i- 

2. Table 2 : 

Experimental data and the averaged values for the ratios i?^', Rj/^. 

3. Table 3 : 

Numerical results for pi, pL + Ap^,, -Rj/^ and Rj/^, — ARj/^ at three values of Ai ; 
Ai, MAX, Ai, MiN and an intermediate value between two extreme values in three scenarios 
(^2 = 2,1,0). 

4. Table 4 : 

Experimental data and the averaged values for the ratios 5", 5**. 
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